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ABSTRACT 

The physics of the "hot spots" on stellar surfaces and the associated variability of accreting magne- 
tized rotating stars is investigated for the first time using fully three-dimensional magnetohydrodynamic 
simulations. The magnetic moment of the star \x is inclined relative to its rotation axis $1 by an angle 
O (we will call this angle the "misalignment angle") while the disk's rotation axis is parallel to f2. A 
| sequence of misalignment angles was investigated, between = 0° and 90°. Typically at small O the 

0^ ■ spots are observed to have the shape of a bow which is curved around the magnetic axis, while at largest 

O the spots have a shape of a bar, crossing the magnetic pole. The physical parameters (density, veloc- 
ity, temperature, matter and energy fluxes, etc.) increase toward the central regions of the spots, thus 
the size of the spots is different at different values of these parameters. At relatively low density and 
temperature, the spots occupy approximately 10 — 20% of the stellar surface, while at the highest values 
of these parameters this area may be less than 1% of the area of the star. The size of the spots increases 
& ■ with the accretion rate. The light curves were calculated for different O and inclination angles of the 

Onj disk i. They show a range of variability patterns, including one maximum-per-period curves (at most of 

angles O and i), and two maximum-per-period curves (at large O and i). At small O, the funnel streams 
may rotate faster/slower than the star, and this may lead to quasi-periodic variability of the star. The 
| results are of interest for understanding the variability and quasi- variability of Classical T Tauri Stars, 

■ millisecond pulsars and cataclysmic variables. 
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1. INTRODUCTION spots. In spite of the large volume of the observational 

T , , . a , , ■ data on variability of these stars, relatively little is known 

In accreting magnetized stars the inflowing matter is •;. ' J 



H . v n j i • 4-1, i- i c iu 4. about the properties of the hot spots 

channelled to regions near the magnetic poles of the star, _ i • , • , m r£ /^,mm^\ i 

r ■ i„ , , „ , , r / , n hoi example, m ciassicai 1 lauri stars (CI IS), hot 

forming hot spots on the stars surface e.g., Ghosh & , ..,.'. . . . ... ... v r ,. 

t u inm tt c cu-u 4. moc n „ ■ j mnn spots radiation is associated with veiling of the contm- 

Lamb 1979; Uchida & Shibata 1985: Camenzmd 1990; r .. . . . . . . ° . r ^.m™ 

Kdnigl 1991; Shu et al. 1994; Hartmann et al. 1994). uum radiation and was analyzed for a number of CTTS 

The hot spots form under conditions in which the star's Herbst k A Koret (1988) and Bouvier et al. (1986, 1993 

,. n i j • i 4. £ i estimated parameters of the not spots assuming a modcf 

magnetic fiefd is strong enough to form a magnetospnere. .,, . . . , . .. , T • „ , , „ 

„ , -i-i i t ,mm • . / TJU4.4-1 with black body radiation, whiic Lamzm (1995), Caivet & 

Examples include classical 1 lauri stars (e.g., Herbst et al. „ „. . „. .' . /„„„„\ v a i-i o 

1986; Bouvier & Bertout 1989; Johns & Basri 1995; 1998; g ullb ™ g i " 8) j ^ullbrmg e \ & \- (2 °° i, i 

Petrnv et al 9001 a h- AWar Rr Ratalha 900?^ rata BaSrl 2000 ) USed models of radiation from a shock wave. 

clvsmic variahles fe p- WickramasWhe Wn Hr Ferrarin The y conclu ded that the area of the star covered with the 

ciysmic vaiiaijics icg., vv icKianiasiiigiic, vvu, c*c rciiaiio . „. n . „ M nnftf . 

inm t- • i n • i mno w innc w nnn\ hot spots m different CI IS varies from 0.3% to 20% de- 

1991; Livio & Prmglc 1992; Warner 1995; Warner 2000), .. ,. ,. , . ,. r , .. 

v i / V,, t, p t u mm rr •• ,i pending on the accretion rate and other factors. In alf 

X-ray pulsars (e.g., Ghosh & Lamb 1979; Irumpcr et al. , , , ,, , , , • , , , , , 

inoc T3-U 4. 4-1 inn7\ j -n- j i / models, however, the hot spots were considered to be no- 

1985; Bildsten et al. 1997) and millisecond pulsars (e.g., ' ,, , '. F , , , 

riu t u i. 4-i nnno\ u • , , i-rp , mogeneous, that is, of the same density and temperature. 

Ghakrabarty, et al. 2003). Ihese objects have different „ °. ' ri ' . . .. . , , , . , , 

- , ,. ' a , , , ,-, i i further, some of the modcis did not take into account the 

dimensions and magnetic fiefd strengths, but the undcrly- „, , ' . . . . »,,... „ , 

r,, , i- ,. jij- j.; effect of imib-darkemng. Additionally, the shape and loca- 

mg physics of tne magnctospneric accretion and not spots . „ , , ° , , J \ , , , 

, • , , . , • -1 tion of the hot spots were not known because they can be 

properties are expected to be simifar. . . . , r ,i t • i i • -i 

m, i j i.i., / . • , r , i i derived onfy from a threc-dimensionaf analysis. Similar 

Ihe observed variability/ quasi- variability of these ob- . , J , ,. . . , J . ,. „ 

i , , r -ii i problems and questions appear during analysis of A-ray 

iects can be due to a number of possibie phenomena m * , , , . .... „, ., . . , , 

li , i i td • nnno n 4- nnno^ pufsars and cataclysmic variabfes. Ihus, it is important 

the magnetosphere (see, e.g., Bouvier 2003; Petrov 2003 . : , , ,, ,. r ,i , , , • , , , ■, 

^ ; r , . \ , . , ... , ' to study the properties of the hot spots m greater dctaii. 

One of the time-scaics may be associated with the hot _ , • r f , , i. <. 

bor analysis of the hot spots we use results of our previ- 
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Fig. 1. — Matter flow near the star and hot spots for the case of a relatively warm disk (T<j = 0.03). The square panels 
show the surfaces of the density in the magnetospheric flows (green color corresponds to p — 0.3 — 0.4, while yellow-green 
color to p — 0.6). Red lines are sample magnetic field lines; the dark-green lines show the sample streamlines of matter 
flow. The circular panels show density distribution in the hot spots for different misalignment angles 0. For = 15°, 30° 
and 45°, the magnetic axis is directed towards the observer, while for = 60°, 75° and 90°, it is directed at an angle 45° 
relative to the observer. It is shown by the crosses in the circular panels. The direction of the rotation axis is shown with 
the open circles. 



ous three-dimensional simulations (Romanova et al. 2003, 
hereafter - R03), and also perform new simulations at a 
lower temperature in the disk and at a variety of parame- 
ters of the star and the accretion flow. 

The main questions which can be answered are: (1) 
What are the shapes of the hot spots? (2) What is the 
distribution of physical parameters (density, temperature) 
in the spots? (3) What is the area covered by the spots 
at different physical parameters? (4) What are the ob- 
served light curves at different misalignment angles and 
different inclination angles of the disk il 

In §2 of the paper we describe the underlying model and 
dimensional examples for CTTS and millisecond pulsars. 
In §3 we discuss the expected physical properties of the 
hot spots. In §4 we calculate the intensity of radiation 
from the hot spots and the associated light curves. In §5 
we discuss dependence of results on different parameters, 
limitations of the model and future work. §6 gives the 
summary of our work. 

2. THE MODEL AND REFERENCE VALUES 

Below we briefly describe the numerical model of R03, 
initial and boundary conditions, and reference values for 
CTTS and millisecond pulsars. 



2.1. The Model 

Disk accretion to a rotating star with an inclined dipole 
magnetic field was investigated with three-dimensional 
MHD simulations. The magnetic moment fi of the dipole 
is inclined to the star's rotation axis by an angle 0. The 
rotation axis of the star coincides with the rotation axis 
of the disk. The value of the magnetic moment p — \fi\ is 
restricted by computer resources. It is difficult to calculate 
in three dimensions large magnetospheres, where the mag- 
netospheric radius R m is much larger than the star's radius 
This is because of the strong variation (oc 1/-R 3 ) of 
the dipole magnetic field. In R03, simulations were done 
for R m /R* w2-5. These values are appropriate for many 
Classical T Tauri Stars (CTTS) and millisecond pulsars. 
For larger magnetospheres, the inner boundary may be 
interpreted as an intermediate layer of magnetosphere. 

A Godunov-type numerical code was used to solve the 
full system of ideal MHD equations in three-dimensions 
space written in a "cubed sphere" coordinate system ro- 
tating with the star (Koldoba et al. 2002; R03). We use 
a reference frame (X, Y, Z) rotating with the star. This 
frame is oriented such that the Z axis is aligned with the 
star's rotation axis and vector fi is in the (X, Z) plane. 
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Fig. 2.— Same as Figure 1 but for cooler disk (Tj = 0.01). 



The quasi-equilibrium initial conditions were proposed 
and tested in axisymmetric model in Romanova et al. 
(2002 - hereafter R02) and in three dimensions in R03. 
These initial conditions take into account initial balance 
of gravitational, pressure gradient and centrifugal forces. 
Although this initial condition does not include the mag- 
netic field, it avoids the possible initial discontinuity of 
the poloidal magnetic field lines at the boundary between 
the disk and corona (corona above the disk rotates with 
the angular velocity of the disk). Test simulations with a 
non-rotating corona show strong magnetic braking and ac- 
cretion of the disk with a speed close to the free-fall speed 
(see Hayashi, Shibata & Matsumoto 1996, some runs of 
Miller and Stone 1997). The new initial conditions de- 
crease this initial magnetic braking dramatically, but not 
completely. There is a residual small magnetic braking 
which determines a slow inward accretion of matter to the 
star. This accretion is used as a source of accreting matter 
in most of simulations. To check the validity of this ap- 
proach, we added simplified "alpha" viscosity to the code, 
where only the largest terms were taken into account (as 
in R02) and performed a set of simulations at = 15° 
and at a variety of a parameters (see §5.1). Simulations 
have shown that the velocity of the accretion flow induced 
by the residual magnetic braking corresponds to viscous 
flow with a « 0.01 — 0.02. These initial conditions allowed 
us to observe for the first time the magnetospheric funnel 
flows from the disk to the star in two and three dimensions 
and to investigate them in detail (see R02, and R03). 

We consider the case when the star rotates relatively 



slowly with angular velocity f2» w O.OAQk*, where ^a'* = 
y/GM/Rl - is the brake -up angular velocity of the star. 
For the case of T Tauri stars, this corresponds to a slow 
rotation rate (T ks 9.4 days for the parameters used in 
R02). This rotation velocity was chosen simply because it 
was the main case of R03. Test cases with faster rotation 
of the star were also performed. 

We consider two main cases: (1) a relatively hot disk 
where the flow to the star is subsonic, and (2) a cooler 
disk, where the flow is supersonic. The important param- 
eters are the initial densities in the disk pd and corona p c , 
and the initial temperatures in the disk Td and corona T c . 
These values are determined at the fiducial point at the 
boundary between the disk and corona (at the inner ra- 
dius of the disk) such that to support pressure balance in 
vertical direction. In two-dimensional simulations (R02) 
we used such parameters that: p c = Q.Olpd, Tg — 0.01T C 
and obtained supersonic funnel streams, because the sound 
speed is relatively low, c s rs O.Ivk- However, in 3D simu- 
lations in R03 we were able to cover only a "softer" range 
of parameters: p c = 0.03pd 7 Td = 0.031^, when the disk 
has higher temperature and c s ks 0.17vk- In this paper we 
use results of R03 and newer runs with lower temperature 
in the disk (c s ~ O.Ivk)- 

The boundary conditions are similar to those in R03. 
At the stellar surface there are "free" boundary conditions 
to the density and pressure. The poloidal components of 
the magnetic field are determined from the fact that the 
star is treated as a perfect conductor rotating at the rate 
fi = Slz. There is however a "free" condition to the az- 
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Fig. 3. — Fraction of the star's surface area / covered by the hot spots with density larger than p for different misalignment 
angles 8. The left panel shows the case of warmer disk, Td = 0.03, while the right panel the case of cooler disk, Td — 0.01. 



imuthal component of the magnetic field: d(RB ( p)/dR = 
so that magnetic field lines have a "freedom" to bend near 
the stellar surface. In the reference frame rotating with 
the star the flow velocity is adjusted such that to be par- 
allel to B at R = i?*. This inner boundary condition 
is valid when the flow is subsonic. In the supersonic case 
there will be standoff shock, but it requires separate inves- 
tigation. From the other side, these boundary conditions 
are valid for any of these cases if the boundary is not the 
surface of the star, but corresponds to some intermediate 
layers of the magnetosphere, while the real surface of the 
star is at smaller radii. Thus, the obtained parameters of 
the hot spots reflect the distribution of these parameters 
in the cross-section of the funnel streams. 

At the outer boundary R = i? max , free boundary condi- 
tions are taken for all variables. We investigated different 
types of boundary conditions for the magnetic field, but 
did not observe a significant difference, because the mag- 
netic field of the dipole decreases very rapidly with dis- 
tance and the difference is much less dramatic compared to 
monopole magnetic field (Ustyugova et al. 1999). We also 
investigated the possible influence of the outer boundary 
conditions by running cases where the simulation regions 
had Rmax/R* = 14, 40, and 194. We find the same results 
except for the case P m ax/P* = 14 where the accretion 
rate decreases too fast because the reservoir of matter in 
the disk is too small. Results for the medium and large 
regions are very close, so that we took R max /R* = 40 as a 
standard size of the simulation region. 

2.2. Dimensionless Variables and Reference Values 

We use the following dimensionless variables: the length 
scale r' = r/i? , the fluid velocity v' = v/v , the density 
p' = p/po, the magnetic field B' = B/B a , the pressure 
p' = p/po, the temperature T' = T/T and time t' = t/to. 
The 'zero' subscript variables are dimensional reference 
values, which are different for different objects. 

The reference values arc determined as follows: the unit 
of distance Rq is taken to be an initial inner radius of 
the disk, R = (Rd)t=o, so that r = 1 at this radius. 
The star has radius P* = 0.35Pq. The reference veloc- 



ity is the Keplerian velocity at Po, = {GM / Rq) 1 / 2 , 
where luq = vq/Rq is the angular velocity. The refer- 
ence time is t = Ro/v - However, in discussing our re- 
sults we measure time in units of the Keplerian period 
of the disk, Po = 27rfo- The reference magnetic field Bo 
is the initial magnetic field strength at r = Rq. The 
reference density is taken to be po = Bq/v 2 . The ref- 
erence pressure is p = pov 2 . The reference tempera- 
ture is To = PofR-p = v 2 /TZ, where TZ is the gas con- 
stant. The reference accretion rate is M = poVqRq- 
The reference energy flux is E = p VqRq. The reference 
value for the black-body temperature of the hot spots is 
(T e ff) = (po^o/c) 1 / 4 , where a is the Stephan-Boltzmann 
constant. 

Subsequently, we drop the primes on the dimensionless 
variables and show dimensionless values in most of the 
figures. Two types of initial conditions for the disk in di- 
mensionless variables are considered: pd = 1, p c = 0.03, 
Td = 0.03 and T c = 1 for the warmer disk, and pd = 1, 
p c = 0.01, T d = 0.01 and T c = 1 for the cooler disk. Below 
we discuss dimensional examples for stars with relatively 
small magnetospheres, CTTS and millisecond pulsars. 

2.3. Reference Values for Classical T Tauri Stars 

Here, we discuss the numerical parameters for a typical 
CTTS. We take the mass of the star to be M* = 0.8M Q , 
and its radius P* = 1.8P Q . The magnetic field at the 
surface of the star is assumed to be _B» = I0 3 G. The 
reference value of length is i?o = (l/0.35)i?» « 2.86i?» w 
3.6 x 10 11 cm. The size of the simulation region corre- 
sponds to i? max w 40i?» w 0.34 AU. The reference ve- 
locity is v ~ 1.9 x I0 7 cm/s. The period of Keplerian 
rotation of the inner radius of the disk is Po ~ 1-38 days. 
We consider that the star rotates relatively slowly, P* = 
9.4 days which corresponds to many CTTSs. Simulations 
for higher rotational velocities did not bring principally 
new results (excluding the propeller stage of course, which 
we do not discuss here). The reference magnetic field 
is P = P*(P»/P ) 3 ~ 42.7 G. The reference density 
is po = 4.9 x I0~ 12 g/cm 3 or n = 3.06 x I0 12 cm~ 3 
which is typical for T Tauri star disks (see, e.g., Hart- 



Fig. 4. — An X — Z slice through the middle of the funnel stream at 8 — 15°. The contour lines show the plane 
cross-section of the density distribution inside the funnel stream. The density changes exponentially from p — 0.2 (blue) 
to p = 2.0 (red). The density in the corona above the disk is p = 0.01 — 0.02. Red lines with arrows show selected magnetic 
field lines. The magnetic moment /x, and the rotation axis fi are shown. 



mann et al. 1998). The reference temperature is To = 
4.5 x 10 6 K which corresponds to typical temperature in 
the corona and temperature Td — 4.5 x 10 4 K in the in- 
nermost part of the disk for the run typical for CTTS 
{Td = 0.01, T c = 1). The reference mass accretion rate is 
M 1.2 x 10 19 g/s k 1.9 x 10~ 7 M Q /year. The refer- 
ence value for the energy flux is Eq w 4.5 x 10 33 ergs/s 
and for black-body temperature: (T e tf)o = 5.0 x 10 3 K. 
The dimensionless values of the intensity in the Figures 7- 
9 are Jw2x 10~ 3 — 1.7 x 10 -2 which corresponds to the 
dimensional values JE w 9.0 x 10 30 — 7.6 x 10 31 ergs/s. 

2.4. Reference Values for Millisecond Pulsars 

We take the mass of the neutron star to be M* = 
1.4 M = 2.8 x 10 33 g, its radius i?„ = 10 km = 10 6 cm, 
and the surface magnetic field B = 10 s G. The reference 
length scale is R w 2.86P* = 2.86 x 10 6 cm. The refer- 
ence velocity is vq = 8.1 x 10 9 cm/s. The reference period 
is Po = 2.2 x 10~ 3 s. In application to millisecond pulsars 
often the frequency is analyzed. The reference frequency 
is Vq = 1/Po ~ 454 Hz. The reference magnetic field is 
B = B*{R t /R ) 3 w 4.3 x 10 6 G. The reference density 
is po — 2.8 x 10 -7 g/cm 3 . The reference temperature is 
3p = 7.8 x 10 11 K. The reference mass accretion rate is 
M w 1.85 x 10 16 g/s w 2.9 x 10~ 10 M /year. The refer- 



ence value for the energy flux is Eq ps 1.2 x 10 36 erg/s. The 
reference black-body temperature is (T e //)o = 4.0 x 10 6 K. 
The dimensionless values of the intensity in Figures 7-9 are 
J w 2 x 10 -3 — 1.7 x 10 -2 which corresponds to dimensional 
values JE„ » 2.4 x 10 33 - 2.0 x 10 34 ergs/s. 

3. PHYSICAL PROPERTIES OF THE HOT SPOTS 

In this section we discuss the connection between the 
shape of the funnel streams with that of the hot spots 
(§3.1), properties of matter along the funnel streams 
(§3.2), and properties of the hot spots (§3.3) 

3.1. Funnel Streams Versus Hot Spots 

The shape and other properties of the hot spots are 
determined by those in the cross-section of the incoming 
funnel stream. In R03 we observed that the shapes of the 
funnel streams are different for different misalignment an- 
gles and may be complicated. At small misalignment 
angles, O <^ 30°, matter typically accretes in two streams. 
For "medium" angles, 30° <, ® <, 60°, the streams are 
often split into several streams. At even larger angles, 
;> 60°, matter again accretes in two streams, but they 
have different shapes compared to those at small 0. The 
streams tend to settle in the particular location, close to 
the fi — Cl plane, after few rotation periods Po- They settle 
faster for cooler disks and large 0. We chose one moment 
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of time t = 5Pq, at which many streams settled, and in- 
vestigate spots in detail for this time. 

Figure 1 shows the funnel streams and corresponding 
hot spots at the surface of the star at time t = 5P a . Only 
a small part of the simulation region is shown in order to 
resolve the inner regions of the funnel streams in greater 
detail. To show the shape of the funnel streams, we have 
chosen for each plot a fixed density level which is slightly 
different for different 6 and varies between p = 0.35 — 0.45 
for smaller O (green color), and p = 0.6 — 0.7 (yellow-green 
color) for larger O. At these density levels the size of the 
spots is relatively large and they occupy approximately 
20% of the star's surface area. One can see that at all 
the density of the spots increases towards their central re- 
gions. These results correspond to simulations of warmer 
disk, T d = 0.03. 

Figure 2 shows similar results, but for the cooler disk, 
Td = 0.01. One can see that the funnel streams and the 
spots are somewhat similar. However, in this case (which 
is taken also for the time t = 5Po) the streams and spots 
are located much closer to their "final" , quasi-equilibrium 
position (which is downstream of the p — fi plane for slowly 
rotating star, - R03), compared to the warmer disk. In 
both cases, of warmer and cooler disks, the spots continue 
to wander around this quasi-equilibrium position. How- 
ever, in the case of the cooler disk the amplitude of these 
oscillations is much smaller. 

We calculated the filling factor of the spots / (the frac- 
tion of the star covered by the spots) for both cases at dif- 
ferent density levels p and different 0. Figure 3 shows that 
at given the value / decreases with p, that is, at larger 
densities the area covered by the spots is smaller. In case 
of warmer disk (left panel), the size of the spots strongly 
depends on 0, specifically at larger densities, while in case 
of cooler disks this dependence is much less prominent. At 
the density levels corresponding to the funnel streams, the 
size of the spots is relatively large and they occupy ap- 
proximately 10 — 20% of the star's surface area. The spots 
are typically smaller for the cooler disk. 

3.2. Parameters Along the Funnel Stream 

To investigate parameters along the funnel stream, we 
consider the case, = 15°, where the pattern of the 
streams is relatively simple. Figure 4 shows the slice 
through the middle of the funnel stream for Td = 0.03. 
One can see that matter accumulates near the magneto- 
sphere forming a relatively dense region - a ring. When 
matter goes to the funnel stream, the density along the 
flow first decreases, but later increases due to the converg- 
ing of the magnetic field lines. Figure 4 shows that the 
density is largest in the interior regions of the stream and 
decreases outward to much lower values. 

Figure 5 shows velocities along «| and across v t the 
stream and Mach number M — v/c s , where v is the total 
velocity. The left panel corresponds to Td — 0.03. One can 
see that in the beginning of the stream, the velocity across 
the stream v t is comparable with v\\, but it decreases while 
approaching the star. The velocity does not reach super- 
sonic values, because the sound speed is relatively high and 
also because the pressure gradient decelerates the flow. In 
case of cooler disk (right panel) the U|| component is larger, 
the Mach number is larger than unity and flow is super- 



sonic. Below we investigate the properties of the hot spots 
for these two cases. 

3.3. Density, Velocity and Matter Flux in the Spots 

Figure 6 shows an example of distribution of different 
parameters in the hot spots for = 30° for cases of sub- 
sonic flow (top two rows) and supersonic flow (bottom two 
rows) . One can see that both, density p and total velocity 
v increase towards the central regions of the spot. The 
distribution of the Mach number M = v/c s is similar to 
that of v. The Mach number reaches values M ^ 2 in the 
middle of the spot in case of cooler disk. 

The matter flux at each point R of the star's surface is 



F M (R) = p n • v 



(1) 



where n = — f is the inward pointing normal vector to the 
star's surface. The total mass accretion rate is 



M 



dfi F M (R) 



where aYl is the solid angle element. Figure 6 shows that 
the Fm distribution is similar to those of density and ve- 
locity. Note, that in both, subsonic and supersonic cases 
the maximum values of Fm are approximately the same, 
which reflects the fact that the accretion rate from the disk 
is approximately the same in both cases. 

3.4. Energy Flux and Temperature in the Spots 

Matter flowing to the star in an accretion stream carries 
both kinetic and thermal energies. The associated energy 
flux through the surface of the star at the point R is: 



F E (R) = p n v ( 2 v2 + w 



(2) 



where, v is the velocity of plasma relative to the surface 
of the star, w — 7(p/p)/(7 — 1) is the specific enthalpy of 
the plasma, and n = — f is the inward pointing normal to 
the surface of the star. 

The detailed physics of the hot spots is complicated. 
To calculate the spectrum of radiation from the hot spots 
an analysis of the radiation transfer is needed (see, e.g., 
Lamzin 1998; Muzerolle, Hartmann, & Calvet 1998; Muze- 
rolle, Calvet, & Hartmann 2001). This paper does not 
include an analysis of the radiation transfer. We obtain 
approximate temperature distributions in the hot spots 
based on overall energy conservation. Namely, we as- 
sume that energy released in the hot spots is due to 
the energy flux of equation (2). At the point R of the 
surface of the star, this energy flux is considered to ra- 
diate as a black body. Thus, Fe(R) = oT 4 ff , where 
a = 5.67 x 10~ 5 (crg/s)/cm 2 /K 4 is a Stcphan-Boltzmann 
constant and T e g is the effective black-body temperature 
of the radiation. Thus, we get 



T cff (R) = 



' p n v 












(7 







"I 1/4 



(3) 



Figure 6 shows the distribution of F^ and T e g in the 
spots. Their shapes are similar to those of the density. 
The distribution of T c g is important for understanding the 
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Fig. 5. — Distribution of velocities along a magnetic field line going approximately through the middle of the funnel 
stream obtained at 6 = 15°. Variable s is the linear distance along this magnetic field line, with s = at the surface of 
the star. Here, W|| is the velocity along the field line, v t is the velocity in the perpendicular direction, M = v/c s is the 
Mach number based on the total velocity of the flow. The left and right panels show the examples of the subsonic and 
supersonic flows correspondingly. 



variability in different spectral bands. The spots are ex- 
pected to be smaller at higher T e s and larger at lower T e g. 
An indication of such distribution was possibly observed 
in the CTTS object BP Tauri, where the area of the hot 
spots was estimated to be different when different meth- 
ods were used. Errico, Lamzin & Vittone (2001) estimated 
the area to be 20% of the surface of the star. They used 
the spectral lines which are thought to originate in the ex- 
ternal regions of the funnel flow. Ardila & Basri (2000) 
estimated the area to be much smaller, 0.3%. However, 
they modelled the UV continuum, which originates in the 
hottest region of the spots. Note that in present analysis, 
we do not take into account the radiation coming from the 
interior of the star. 

4. THE INTENSITY OF RADIATION AND LIGHT CURVES 

For calculation of the radiation intensity from the hot 
spots, we suppose that all of the kinetic energy flux Fg(R 
goes into radiation. The energy, radiated by a unit area 
of the spot is 

F E (R)= J dn/(R,m), (4) 

cos0>O 

where /(R, m) is the intensity of the radiation from a unit 
area into the solid angle element dfl in the direction m 
with m • f = cos 9. The condition cos 9 > corresponds to 
the upper half space above the star's surface. For speci- 
ficity we consider /(R, m) = A(R)cos9. Thus we obtain 
the luminosity of unit area at point R, 

tt/2 

Fe(R) = 2irA(R) J sin 9 cos 9d9 = ttA(R) , 
o 

and 

/(R,m) = -F E (R)cos<9 . (5) 

Consider now the intensity of radiation seen by a distant 
observer. The unit vector from the star to the observer is 



denoted k. This vector makes an angle (3 with the normal 
to the star's surface, that is, cos/3 — k • f . The intensity 
of radiation in the direction k is obtained by integrating 
over the stellar surface, 

J = J dS /(R, k) = i J dS Fe(R) cos /? , (6) 

cos,3>0 

where dS is an element of the surface area of the star. The 
condition cos /3 > or (3 < tt/2 corresponds to only the 
near side of the star being visible. The rotation of a star 
with hot spots leads to variations in the observed emission 
J. 

4.1. Variability Curves at Fixed position of Spots 

Hot spots constantly change their shape and position. 
However, in many situations the changes are small (e.g., in 
case of the cooler disk), so that as a first step we suppose 
that the spots have fixed location at the surface of the 
star and variability is connected with rotation of the star 
relative to the fixed observer. This approach helps to un- 
derstand different features of the light curves, associated 
with structure and location of the hot spots at different 6. 
Examples of actual variability curves are shown in §4.2. 

We consider one of the initial conditions (warmer disk, 
Td = 0.03), take hot spots at a single moment of time, 
t = 5Po (as in Figure 1), and calculate from equa- 
tion (6) the observed intensity J(t) as the star rotates. 
Such light curves were calculated for misalignment angles 
9 = 15°, 45°, and 6 = 75°, and for different inclination 
angles of the disk i with respect to the line of sight. The 
inclination angle of the disk is the angle between the di- 
rection of fi and the direction to the observer k. That is, 
cosi = fi • r . The results are shown at Figures 7-9. 

Figure 7 shows results for a relatively small misalign- 
ment angle = 15°. The top panel shows the light curves 
for different inclination angles i. The bottom panel shows 
the orientation of the star during one rotation period. The 
phases are shown for two cases: at i = 30°, which is a 
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Fig. 6. — Distribution of different physical parameters in the hot spots for a misalignment angle = 30° and for two 
types of initial conditions in the disk and corona, p is the density, v is the total velocity, T e ff is the effective temperature, 
M is the Mach number of the flow, Fm and Fe are the fluxes of mass and energy. The top two rows show the case of 
higher temperature in the disk and subsonic flow, while the lower two rows show the case of lower temperature in the disk 
and supersonic flow. 



probable inclination angle, and at i = 90° which is edge 
on. 

Of course for i = 0° there is no variability because the 
observer sees only one hot spot which rotates around the 

axis. For larger inclination angles, i — 15°, i — 30°, 
i = 45°, dips appear in the light curves which increase 
with increasing of i. These dips are connected with the 
fact that the hot spots come closer to the edge of the star 
and are partially obscured. Under these conditions the 
star will be observed as a variable star with one maximum 



per period. For i — 60°, the light curve becomes asym- 
metric and this asymmetry increases for larger inclination 
angles i. This is connected with the fact that the second 
spot, which was invisible at smaller i, appears and starts to 
contribute to the luminosity. Thus, two intensity maxima 
per stellar period are observed. The two maxima become 
of equal amplitude at i = 90°. 

Figure 8 shows similar plots, but for the misalignment 
angle 9 = 45°. The light curves in this case are similar 
to those at smaller 0. However, the shorter time-scale 
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Fig. 7. — Light curves calculated for small misalignment angle O = 15° and for different inclination angles i. Time t is 
measured in periods of rotation of the star P. The two bottom pictures show the distribution of energy flux F(R) in the 
hot spots at different phases ip of the star seen by observer at inclination angles i = 30° and i = 90°. Phases are shown 
for one period of rotation. The magnetic and rotation axes are marked with a circle and a cross correspondingly. 



variability appears at large inclination angles at i > 75°. 

Figure 9 shows the case of a high misalignment angle, 
9 = 75°. For this and larger values of 0, the double max- 
ima (half-period) variability appears for many inclination 
angles starting from i J> 30°. 

The figures 7-9 show intensity in dimensionlcss units. 
This value in dimensional units will be JEq, where Eq is 
the dimensional energy flux (see §2.3 and §2.4). The value 
JEq gives the energy flux from the spots radiated to the 



direction of the observer. This value will determine the ob- 
served flux, when the distance to the object will be taken 
into account. Thus, the X-ray flux observed e.g. from mil- 
lisecond pulsars will be proportional to J and the expected 
variability curves will be similar to those at the Figures 7- 
9. However, for CTTS it is common to show the light 
curves on a logarithmic scale, using the standard stellar 
magnitude values. As long as this analysis can be appli- 
cable to stars at different distances, we accept the random 
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Fig. 8. — Light curves calculated for medium misalignment angle 
features of the plot are the same as in Figure 7. 



45° and for different inclination angles i. All 



null-point for calculation of m which will be different for 
different CTTS. We suppose that intensity J = 0.1 cor- 
responds to m = 12, and calculate the light curve using 
the standard formulae: m = 12 — 2.512(1. + logj). Figure 
10 shows light-curves in stellar magnitudes for different 9 
and i. Note that the curves have different shapes which 
may be used for estimating values of and i. 



4.2. Variability Curves for Changing Spots 



Here, we present samples of actual variability curves. 
We take the data for all moments of time from the numer- 
ical simulations. The star rotates slowly, with f2* = 0.19, 
so that one rotation of the star corresponds to 5.26 rota- 
tions of the inner radius of the disk P» « 5.3p). Thus, 
to show the variability pattern during few rotations of the 
star, we need pretty long runs. Here, we show sample re- 
sults from several relatively long runs with t = 15 — 16Po. 
Figure 11 shows sample light curves for the case of the 
cooler disk, and O = 15° and = 45°. One can see 
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Fig. 9. — Light curves calculated for a large misalignment angle O = 75° and for different inclination angles i. All features 
of the plot are the same as in Figure 7. 



that the "real" light curves have similar features to the 
curves for fixed hot spots. However, the variability pat- 
tern departs from the exact variability patterns shown at 
Figures 7 and 8. This is because the hot spots continue 
to change their shape and position. Figure 12 shows vari- 
ability curves for larger misalignment angles, = 75° and 
Q = 90° (for warmer disk). One can see that at these 
misalignment angles, the light curves are closer to those of 
Figure 9, because the spots have almost fixed positions on 
the stellar surface. 



5. DISCUSSION 

In this paper we fixed parameters of the star and the 
flow. Below we discuss dependence of the results on fi* 
and on the accretion rate M. In §5.2 we discuss the limi- 
tations of the model and future work. 

5.1. Dependence of Results on Parameters of the Star 
and the Disk 

Q-Dependence: For the results given here we have taken 
a relatively low angular velocity, f2» = 0.19. This corre- 
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sponds to CTTS with a period P = 9.4 days. We did 
additional simulations for an intermediate angular veloc- 
ity, 0* = 0.35 (P = 3.3 days), and a rapidly rotating 
star, fi* = 1.0 (P = 1.8 days). We observed that the hot 
spot shapes are similar for slowly and for rapidly rotat- 
ing stars. For the intermediate angular velocity, the bow 
shape is less prominent. The bow shape of the hot spots 
reflects the typical shape of the funnel streams which have 
a thickness that is small compared with their width. This 
is typical for relatively small inclination angles, G <J 30°, 
where the stream should "climb" an appreciable distance 
above the equatorial plane before accreting to the vicinity 



of the magnetic pole. At small inclination angles, the bow 
shape is natural, because it may be considered as a part 
of the narrow cylindrical ring which occurs for O = 0. 

M -Dependence: We studied the dependence of the size 
of the spots on the accretion rate M. We introduced an 
alpha-type viscosity analogous to one used in R02 and per- 
formed a set of runs at different values of accretion rate M 
(different a). Figure 13 shows that the fraction / of the 
star covered by the hot spots increases with the accretion 
rate for a range of density levels, p = 0.4, 1, 2. Thus, the 
size of the spots increases as the accretion rate increases. 
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Fig. 12. — Light curves calculated for actual time-dependent hot spots for misalignment angles 9 = 75° and 9 = 90°, 
and for different inclination angles i, for the case of the warmer disk, T d = 0.03. 



This is in accord with recent results of Ardila and Basri 
(2000), who analyzed the UV variability of the CTTS BP 
Tau and have shown significant correlation between the ac- 
cretion rate and the filling factor of the shocks. In R03 we 
noticed that at a larger accretion rate, the streams become 
wider and cover a larger area of the magnctosphcrc. This 
may possibly explain the observed variation of the shapes 
of magnctospheric lines with the accretion rate (Muzerollc 
et al. 2001). 

5.2. Limitations of the Model and Future Work 

Variability in different spectral lines in the accreting 
magnetized stars may be associated with different regions: 
inner regions of the disk, magnetospheric streams, hot 
spots, or outflows. In this paper we analyze only the vari- 



ability due to the rotation of the hot spots, while the vari- 
ability associated with other regions will be investigated 
in the future work. The light from the hot spots may be 
also obscured by magnetospheric streams or by the warped 
inner regions of the disk if the star is approximately edge- 
on (e.g., Bouvicr, ct al. 1999, 2003). This paper does not 
include the effect of obscuration. The warping of the inner 
regions of the disk (Aly 1980, Lipunov & Shakura 1980, 
Lai 1999) was not observed in the simulations, but a spe- 
cial set of simulations will be done for investigation of this 
possible phenomenon. 

The magnetic field of the CTTS or millisecond pulsars 
may not be a pure dipole field (Safier 1998; Smirnov et al. 
2003). This can lead to more complicated geometry of the 
hot spots and more complicated variability patterns. How- 
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Fig. 13.— The fraction of the star covered by the hot spots at the density levels: p = 0.4, p = 1, and p = 2 as a function 
of the accretion rate M. Squares show separate runs at different a parameters of the disk viscosity. 



ever, the detailed analysis of the photometric and spectral 
variability of a number of CCTS have shown that most of 
the observed features can be explained by models with a 
dipole magnetic field (e.g., Muzerolle, Hartmann, & Cal- 
vet 1998; Petrov, et al. 2001a, b; Alencar, Johns-Krull, & 
Basri 2001; Bouvier et al. 2003). This paper considers 
only the case where the star's field is a dipole field. Non- 
dipolar field geometries will be investigated separately. 

We observed that the hot spots may rotate more rapidly 
or more slowly than the star, which will lead to quasi- 
periodic oscillations in the light curves and may explain 
QPOs observed in CTTS (see Smith, Bonnell, & Lewis 
1995) and millisecond pulsars (see, e.g., Chakrabarty et 
al. 2003). QPO variability in the disk-magnetized star 
systems may also be associated with oscillations of the in- 
ner radius of the disk (e.g., Goodson, Winglee & B6hm 
1997). We plan to further investigate the quasi- variability 
associated with such phenomena. 

6. SUMMARY 

Disk accretion to a rotating star with a misaligned 
dipole magnetic field has been studied further by three- 
dimensional MHD simulations. This work focuses on the 
nature of the "hot spots" formed on the stellar surface due 
to the impact of two or more funnel streams. We investi- 



gated the shape and intensity of the hot spots for different 
misalignment angles between the star's rotation axis fi 
and its magnetic moment (i. Further, we calculated the 
light curves due to rotation of the hot spots for different 
angles i between the observer's line-of-sight and fi. The 
main results are the following: 

1. For small inclination angles, < 30°, the hot spots 
typically have a shape of a bow which is bent around the 
magnetic pole. At large inclination angles, ;> 60°, the 
shape becomes bar-like. Often a spot on a given hemi- 
sphere splits to form two spots, which reflects the split- 
ting of the funnel stream into two streams. The secondary 
stream is typically weaker than the main stream so that 
one spot is much larger than the other. 

2. The density, temperature, matter flux and other pa- 
rameters increase towards the central regions of the spots, 
so that the spots are larger at lower temperature/density 
and smaller at larger temperature/density. They cover 
about 10 — 20% of the area of the star at the density level 
typical for the external regions of the funnel streams (see 
Figures 1 and 2). The size of the hot spots increases with 
the accretion rate. 

3. The spots have tendency to be located close to the p — Ci 
plane. They tend to be located downstream of this plane 
if a star rotates slowly (i.e., the inner region of the disk 
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and the foot-points of the stream rotate somewhat faster 
than the star), or upstream, if a star rotates relatively 
fast. The spots wander around their "favorite" position. 
The amplitude of wandering is smaller in case of the cooler 
disk. 

4. The calculated light curves reveal the following fea- 
tures: 

(a) . The light curve has one peak per one period of rota- 
tion of the star and the shape is approximately sinusoidal. 
This is typical for small and medium misalignment angles, 
9 <; 45°, and inclination angles i £ 60°. (b). The light 
curve has two peaks per period of rotation. This is typical 
for all 9 if inclination angle is large, i > 75°. At very large 
misalignment angles, 9 ^ 60° — 70°, the double-peak curve 
is typical for wide range of inclination angles, i ;> 30°. 

5. The variation of the shape and location of the spots will 
lead to departure from the exact variability and to quasi- 
variability. At small misalignment angles, 9 < 30°, the 



streams (and hot spots) may rotate with velocity different 
from that of the star (R03) thus leading to quasi-periodic 
oscillations. 
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